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ABSTRACT 

We analyze RXTE and Swift observations of SWIFT J1729. 9-3437 after 
its outburst from 2010 July 20 to 2010 August 12. We calculate a spin fre- 
quency and spin frequency derivative of 1.8734(8) x 10 -3 Hz and 6.42(6) x 
10 -12 Hz/s respectively from the quadratic fit of pulse arrival times. The 
quadratic fit residuals fit well to a circular orbital model with a period of 
15.3(2) days and a mass function of about 1.3M©, but they can also be ex- 
plained by a torque noise strength of 6.8 x 10~ 18 Hz sec~ 2 . Pulse profiles 
switches from double-peaked to single-peaked as the source flux continues 
to decrease. We find that the pulse shape generally shows no strong energy 
dependence. The hardness ratios reveal that the source becomes softer with 
decreasing flux. We construct a single spectrum from all the available RXTE 
and Swift observations. We find that adding an Fe line complex feature around 
6.51 keV slightly improves the spectral fit and that this feature is more likely 
to originate from the source rather than the galactic ridge. From the pulse 
phase spectral analysis, it is shown that that photon index and folding energy 
of the high energy cut-off vary with varying pulse phase. 

Key words: X-rays: binaries, pulsars: individual: SWIFT J1729.9— 3437 , 
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1 INTRODUCTION 

SWIFT J1729.9— 3437 is a transient X-ray pulsar discovered during all-sky monitoring of 
Swift Burst Alert Telescope (BAT) on 2010 July 13 (Markwardt et al. 2010a). At the same 
time, Rossi X-ray Timing Explorer (RXTE) Proportional Counter Array (PCA) monitoring 
scans of the galactic center region confirmed a gradual increase in flux from the position 
consistent with the Swift position of the source (RA = 262.53746, DEC = —34.61239). 
Consecutive RXTE-PCA pointings identified SWIFT J1729.9-3437 as a pulsar with ~530 
s pulsations (Markwardt et al. 2010b). Markwardt et al. (2010b) also suggested that the X- 
ray spectrum of the source is compatible with basic X-ray pulsar spectra, modeled by an 
absorbed power law with a high energy cut-off. 

In this paper, we analyze RXTE and Swift observations of SWIFT J1729.9— 3437 and 
examine the spectral and timing properties of its outburst from 2010 July 20 to 2010 August 
12 (MJD 55397 — MJD 55420). In Section 2, we mention the observations that we analyze. 
In Section 3 and 4, we present our timing and spectral analysis results. In Section 5, we 
discuss our results and conclude. 

2 OBSERVATIONS 
2.1 RXTE 

11 pointing RXTE-PCA observations of SWIFT J1729.9-3437 are between 2010 July 20 
and 2010 August 8. The total exposure time is about 42 ks. PCA consisted of five co-aligned 
identical proportional counter units (PCUs) (Jahoda et al. 1996) each having an effective 
area of approximately 1300 cm 2 . The energy range of PCA was from 2 to 60 keV, with 
an energy resolution of 18 per cent at 6 keV. The field of view (FOV) at full width at 
half- maximum (FWHM) was about 1°. Generally, various PCUs turned off during PCA 
observations. The response of PCU0 and PCU1 were not well known due to loss of their 
propane layers. Furthermore, the response of top anode layers was modeled better than 
the other layers. Although the number of active PCUs during the observations of SWIFT 
J1729.9— 3437 varied between one and three, we select only data from the top anode layer 
of PCU2 because of the reasons mentioned. 

The standard software tools of HEAS0FT v. 6. 11 are used for the PCA data analysis. 
Data is filtered by selecting times when the elevation angle is greater than 10°, the offset 
from the source is less than 0.02° and the electron contamination of PCU2 is less than 
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0.1. The background extraction for spectra and light curves are done by using the latest 
PCA background estimator models supplied by the RXTE Guest Observer Facility (GOF), 
EPOCH 5C. During the extraction of spectra Standard2f mode data is considered with 128 
energy channels and 16-s time resolution. Relevant response matrices for spectra are created 
by PCARSP V . 1 1 . 7 . 1 . Furthermore we construct pulse phase resolved spectra with the tool 
FASEBIN Rev . 1 . 8 by using Good Xenon mode event data and 125 [is time resolution event 
files (E_125us_64M_0_ls) when Good Xenon mode data is not available. From these spectra, 
we generate energy resolved pulse profiles by obtaining the count rates per phase bin with 
the tool FBSSUM. The Is binned 3-20 keV light curves are produced from Good Xenon and 
E_125us_64M_0_ls events. 

2.2 SWIFT 

After the discovery of SWIFT J1729.9-3437 by Swift-BAT (Barthelmy et al. 2005) all-sky 
monitoring on 2010 July 13 (Markwardt et al. 2010a), a total of 11 follow-up X-ray Telescope 
(XRT; Burrows et al. 2005) observations (each ~3 ks) were carried out between 2010 July 
20 and 2010 August 12. XRT is a focusing instrument on board the Swift satellite (Gehrels 
et al. 2004) which operates between 0.2—10 keV energy range. CCD— 22 detector of XRT 
has an effective area of 110 cm 2 at 1.5 keV with a spatial resolution of 18" and the FOV of 
23.6' x 23.6'. Operation mode of XRT switches between photon-counting (PC), imaging and 
timing; depending on the brightness of the observed source. Pointing observations of SWIFT 
J1729.9— 3437 are in the PC mode. Screened event files are produced by the XRTDAS v . 2 . 7 . 
standard data pipeline package XRTPIPELINE v. 0.12.6. Standard grade filtering (0-12) is 
applied to PC mode data. 

The spectral extraction is carried out by filtering appropriate regions using XSELECT v . 2 . 4b. 
For the observations in which the XRT count rate is high enough (above 0.5 cts/s) to pro- 
duce event pile-up, source region is selected to be annular. We compared the observed and 
nominal point spread function (PSF) (Vaughan et al. 2006) and decided the size of the core 
affected in order to determine the radius of the inner circle of the annulus. For the brightest 
observations an exclusion region of radius ~ 9" is sufficient to correct pile-up. For low count 
rate observations a circular region of radius 47" is selected for the source spectral extrac- 
tion. Source regions are centered on the position determined with XRTCENTROID v. 0.2. 9. 
A circular source-free region with 140" radius is selected for background spectral extrac- 
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Figure 1. 530 s binned light curves from Swift-XRT (0.3-10 keV; upper panel) and RXTE-PCA (3-20 keV, PCU2 top 
layer; lower panel) observations. 

tion. Resulting spectral files are grouped to require at least 50 counts per bin using the 
ftool GRPPHA v.3.0.1, for the x 2 statistics to be applicable. We used the latest response 
matrix file (version v013) and created individual ancillary response files using the tool 
XRTMKARF v. 0.5. 9 with the exposure map produced by XRTEXP0MAP v. 0.2. 7. Spectral 
analysis is performed using XSPEC v . 12 . 7 . 0. 

For the timing analysis, arrival times of XRT events are first corrected to the solar 
system barycenter by using the tool BARYC0RR v. 1 . 11. Then, background-subtracted XRT 
light curves in the 0.2—10 keV energy range; corrected for pile-up, PSF losses and vignetting; 
are extracted by using the highest possible time resolution (2.51 s) for PC mode data. 
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3 TIMING ANALYSIS 
3.1 Timing Solution 

For timing analysis, we use Is binned RXTE-PCA and 2.5s binned SWIFT-XRT light curves 
of the source. To give an idea of the temporal coverage, in Fig. HJ we present 530 s binned 
versions of these light curves. These background subtracted light curves are also corrected 
to the barycenter of the solar system. 

In order to estimate pulse periods of the source, RXTE-PCA time series is folded on 
statistically independent trial periods (Leahy et al. 1983). Template pulse profiles are formed 
from these observations by folding the data on the period giving maximum x 2 ■ Then the 
barycentered Swift time series are also folded over the best period obtained from RXTE. 
Pulse profiles consisting of 20 phase bins are represented by their Fourier harmonics (Deeter 

© 2013 RAS, MNRAS OOP. HUTU 



SWIFT J 1729. 9- 34 3 7 5 




Figure 2. Pulse phases and their residuals after a quadratic fit. Crosses and triangles represent data from RXTE-PCA and 
SWIFT-XKT respectively. Solid line in the bottom panel corresponds to the circular orbital model with the parameters listed 
in Table Q] 

& Boynton 1985). Using cross-correlation of pulses between template and in each observation, 
we obtain pulse arrival times. 

We are able to connect all pulse arrival times of the observations in phase over the whole 
observation time span. We find that it is possible to fit the pulse arrival times to the second 
order polynomial, 

50 = <50 o + 6u(t - t ) + -u(t - t f (1) 

where 5<p is the pulse phase offset found from the pulse timing analysis, t a is the epoch for 
folding; 8(j> is the residual phase offset at t Q ; 8v is the correction to the pulse frequency at 
time to; v, being pulse frequency derivative, is the second derivative of the pulse phase. In 
Figure El we present pulse phases and residuals of this quadratic fit. From the quadratic fit, 
we obtain the timing solution which is presented in Table [TJ 

We find that these residuals fit well to a sinusoidal function with a period of 15.3 ± 0.2 
days (see bottom panel of Figure 12]) . This corresponds to a circular orbital model, parameters 
of which are listed in Table [TJ 

Quadratic residuals can alternatively express the noise process due to random torque fluc- 
tuations (Bildsten et al. 1997, Baykal et al. 2007). In order to estimate the noise strength, we 
fit cubic polynomial to residuals of pulse arrival times. The observed time series is simulated 
by the Monte Carlo technique for a unit white noise strength defined as Po{f) = 1 and 
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Table 1. Timing Solution and Orbital Model of SWIFT J1729. 9-3437. 



Timing Solution Parameter 


Value 


Timing Epoch (MJD) 


55393.50(6) 


Frequency (Hz) 


1.8734(8) x 10" 3 


Frequency Derivative (Hz.s -1 ) 


6.42(6) x 10" 12 


Circular Orbital Model Parameter 


Value 


Orbital Epoch (MJD) 


55395.700(6) 


^ sini (lt.s) 


65(3) 


Orbital Period (days) 


15.3(2) 
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Figure 3. 26-s binned light curves of three different RX TE— PCA observations (3—20 keV, PCU2 top layer). Time values are 
converted to phases (or time/pulse period) for arbitrary observation epochs. The IDs and mid-times of the observations are A) 
95044-05-02-02 at MJD 55404.6, B) 95044-05-02-03 at MJD 55406.7, C) 95044-05-03-00 at MJD 55408.5. 

fitted to a cubic polynomial in time. Then the square of the third order polynomial term 
is divided into the value from Monte Carlo simulations (Deeter 1984, Cordes 1980). Torque 
noise strength is obtained as 6.8 x 10 -18 Hz sec -2 . This value of noise strength estimate is 
comparable with those of other accretion powered sources such as wind accretors e.g. Vela 
X-l, 4U 1538-52 and GX 301-2 which has the values changing between 10" 20 - 10" 18 Hz 
s~ 2 (Bildsten et al. 1997). Her X-l and 4U 1626-67, which are disc accretors with low mass 
companions, have shown pulse frequency derivatives consistent with noise strengths 10~ 21 
to 10 -18 Hz sec -2 (Bildsten et al. 1997). Therefore resiudals of quadratic fit can also be 
associated with torque noise fluctuation of the source. 
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Figure 4. Energy resolved pulse profiles of RXTE-PCA observations on MJD 55400.78 (left; ID:95044-05-02-00) and MJD 
55408.46 (right; ID:95044-05-03-00). The unit of y-axis of each plot is normalized count from the specified energy interval. The 
dashed line represents mid-point of the secondary peak. 



3.2 Pulse Profiles and Hardness Ratios 

A double peaked pulse profile is observed in the first five RXTE— PCA observations of 
SWIFT J1729.9— 3437 (see panel (A) of Fig. |3]); however one peak loses its intensity starting 
from the middle of the observation on MJD 55406.7 (see panel (B) Fig. [3]) as the source 
flux continues to decrease after the burst. The peak value of 2—10 keV unabsorbed flux is 
3.04 x 10~ 10 ergcm _2 s _1 on MJD 55398.7, during its gradual decrease it reaches to 1.96 x 
iQ-io er g cm -2 g -i on mjj} 55408.5, when the shape of the pulse profile changes (see panel 
(C) Fig. [3]). The minimum flux observed for the last RXTE— PCA observation on MJD 
55416.4 is 1.36 x 10~ 10 ergcm -2 s -1 . These flux values are calculated from the model flux of 
individual spectral fitting of each observation. 

To search for a possible energy dependence of the pulse shape change, we construct pulse 
profiles in five energy bands, i.e. 3—8 keV, 8—13 keV, 13—18 keV, 18—25 keV and 25—60 keV. 
Two examples for energy resolved pulse profiles, one for double peaked and one for single 
peaked, are given in Fig. HJ We find that 8—13 keV and 13—18 keV pulses are stronger (they 
have higher pulse fraction) than the 3—8 keV and 18—25 keV pulses in all observations. 
In 25—60 keV the pulse fraction noticeably drops. The pulse shape shows no strong energy 
dependence except for the two observations on MJD 55400.78 (ID:95044-05-02-00 see Fig. Hj) 
and MJD 55404.63 (ID:95044-05-02-02). In these two exceptions the secondary peak around 
pulse phase 0.3 loses its intensity at the 18—25 keV energy band. 
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Figure 5. Hardness ratios of energy resolved light curves. Daily averaged time evolutions are plotted on the left, whereas 530 
sec binned hardness ratios over 3-25 keV PCU2 count rates are plotted on the right. Specified abbreviations on y-axis of each 
panel stand for HR1: 8-13 keV/3-8 keV, HR2: 18-25 keV/13-18 keV. 

We construct hardness ratio plots from the energy resolved light curves. Daily averaged 
count rates from 3—8 keV, 8—13 keV, 13—18 keV and 18—25 keV light curves are used to 
plot hardness ratios HR1: 8-13 keV/3-8 keV and HR2: 18-25 keV/13-18 keV over time 
(see left panels of Fig. [5]). HR1 shows a noticeable evolution with respect to time, it remains 
almost constant at a value ~ 1.1 until MJD 55406.7 and a gradual decrease starts after the 
sixth observation. The interval with constant HR1 correspond to times when the pulse profile 
is double peaked, where as the decreasing interval of HR1 coincides with the times when 
the pulse profile is single peaked. Nevertheless, it should be noted that the pulse profiles 
show no significant variation in the corresponding energy bands (see Fig. H]). Furthermore, 
530 s binned hardness ratios are plotted over the total count rate in 3—25 keV band (see 
right panels of Fig. [5]). A noticable correlation is again observed for HR1. As the source flux 
decreases during the ongoing decline of the outburst, the emitted radiation becomes softer. 



4 SPECTRAL ANALYSIS 
4.1 Overall Spectrum 

A preliminary spectral analysis for the first pointed observations of SWIFT J1729.9— 3437 
was performed by Markwardt et al. (2010b). In this paper, we extend the spectral study by 
using all the available observations of the source defined in Section 2. Basically, the spectra 
can be modeled by a power law with a high energy cut-off and photoelectric absorption 
as it is suggested. Among the several models that describe the cut-off power law, the best 
fit is achieved by CUT0FFPL model in XSPEC. An additional Gaussian component is also 
required for a weak Fe emission line around 6.4 keV. During the simultaneous fitting of 
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Figure 6. Simultaneous spectral fitting of Swift- 
data, and its best fit with WABSx (CUTOFFPL+GAU) 
the lower panel. 
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(X 2 = 1-03 



Energy (keV) 

9.3 keV; grey) and RXTE-PCA (5-25 keV; black) data. The 
solid line) are shown in upper panel, the residuals are given in 



Table 2. Best fit ( WABSx (CUTOFFPL+GAU)) spectral parameters for the simultaneous fitting of Swift-XRT and RXTE-PCA 
data shown in Fig[6] Ci and C2 are the multiplicative constant factors for different instruments. All uncertainties are calculated 
at the 90 per cent confidence level. 



Parameter 






Value 


Ci 


for XRT 




1.00 (fixed) 


c 2 


for PCA 




1 ccr + 0.04 
i - oo -0.04 


n H 


(10 22 cm" 2 ) 




8 27+ - 37 
° ' -0.36 


r 






, 9 o+0.07 


r Norm. 


(10 -2 cts cm~ 2 s~ 




2A2±°f B 


Efold 


(keV) 




16.50ti : ™ 


Ef e 


(keV) 




D - ol -0.03 


"Fe 


(keV) 






Fe Norm. 


(10~ 4 cts cm -2 s - 


X ) 


4.68tg;« 


Unabs. F . 3-9.3 kcV 


(10 — 10 erg cm -2 s 




9 ns+ 017 

Z - UO -0.024 


Unabs. F 5 _ 25 kcV 


(10 -10 erg cm -2 s 


- 1 ) 


nf ,+0.004 
°- uu -0.016 



Reduced x 2 (dof) 1.09 (187) 



Swift-XRT and RXTE-PCA spectra we included a multiplicative constant factor in the 
model to account for the normalization uncertainty between the two instruments. The data 
and its best fit are plotted in Fig. [6] and the corresponding spectral parameters are given in 
Table El 

The energy ranges for the simultaneous spectral analysis are initially selected to be 
0.3—9.3 keV for XRT spectrum and 3—25 keV for PCA spectrum. Individually these spectra 
have similar shapes that can be modeled with same models, apart from an offset in absolute 
flux calibration. However individual modeling of the data from different instruments yield 

© 2013 RAS, MNRAS OOO.HTfTH 



10 §. §ahiner, S. Q. Inam, M. M. Serim, and A. Baykal 

different absorption parameters due to different instrumental band-passes (Markwardt et al. 
2010b). During the simultaneous fitting trials we observe large residuals for the first energy 
bins of the PCA spectrum although the fit is adequate for XRT spectrum. Therefore we 
exclude energies below 5 keV for PCA spectrum, since XRT has more spectral energy bins 
in soft X-rays. 

The FOV of PCA is large and SWIFT J1729.9-3437 is near the Galactic ridge. The 
thermal emission from the ridge is usually known to contaminate the spectral count rates 
of the source. During the analysis of the first PCA observation, the weak line at ~6.6 keV 
reported by Markwardt et al. (2010b) was found suspicious to be a contamination, since 
it could not be resolved in the spectrum of the first XRT observation. We tried to handle 
this issue in overall PCA spectrum with fixed additive Galactic ridge models. We confirm 
that the flux of the source is two orders of magnitude bigger than the flux of the Galactic 
ridge (Valinia & Marshall 1998) therefore we conclude that the contamination could be small 
enough to be neglected. Furthermore; addition of the Gaussian model component for the 
Fe line improves the individual fit of the 0.3—9.3 keV XRT spectrum after combining data 
from all XRT observations by reducing the reduced x 2 from 1.08 (x 2 /d.o.f. = 159.0/147) to 
0.98 (x 2 /d.o.f. = 143.8/146). The F-test probability that this improvement is achieved just 
by chance is 1.3 x 10~ 4 . Therefore we suggest that the Fe emission originates from SWIFT 
J1729.9— 3437. The line energy found for the source is shifted from the neutral value of Fe 
Ka (~6.4 keV), which may be a consequence of an excess of ionization. 

In accretion powered pulsars, the Fe Ka line is produced by the reprocessing of the hard 
X-ray emission in the relatively low ionized and cool matter surrounding the pulsar. The 
correlation of the iron line intensity and the continuum intensity in many sources (e.g. Vela 
X-l: Ohashi et al 1984, GX 301-2: Leahy et al. 1989) is specified to be the evidence for the 
fluorescence of X-ray continuum by the cold material. As the continuum intensity increases 
the illumination of the matter increases and the line intensity strengthens. In the individual 
analysis of each RXTE observation of SWIFT J1729.9— 3437; the normalization parameter 
of the Gaussian fitted for Fe emission shows a gradual decrease from 8.5 x 10 -4 cts cm" 2 s _1 
to 5.2 x 10~ 4 cts cm -2 s _1 in correlation with the source flux. This correlated change also 
confirms the origin of the Fe emission. 
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4.2 Pulse Phase Resolved Spectra 

We construct pulse phase resolved spectra of SWIFT J1729.9— 3437 from the first three 
RXTE— PCA observations. There are two motivations to the data selection. First, the bright- 
est of the observations are selected to ensure the best signal-to-noise ratio possible. Second, 
the selected observations are the only ones that have 256 channels Good Xenon event mode 
data available, which is required for the FASEBIN tool. The timing solution found for the 
source in Section 3.1 is appended to the timing files of FASEBIN tool for correct phase- 
binning. The pulse period is divided into five equal segments resulting in ~2.1 ks exposure 
time for each phase bin spectrum. 

We model 3—25 keV spectra with the same models used for the overall spectrum in the 
previous section. The Hydrogen column density and Fe line energy are fixed during 
the analysis. The photon index and folding energy of the high energy cut-off are found to be 
varying within pulse phases (see Fig. [7j). The phase dependence of photon index is in anti- 
correlation with the pulsed flux. The folding energy shows a similar trend as the photon 
index, which means softer spectra have higher cut-off energies. Our results indicate that 
the softest emission is observed at the un-pulsed phase around 0.05, where as the pulsed 
emission is harder. The interpulse phase around 0.65 is softer than the phases of pulse peaks. 

5 SUMMARY AND DISCUSSION 

In this paper, we study timing and spectral properties of SWIFT J1729.9— 3437 using RXTE 
and Swift observations of its outburst from 2010 July 20 through 2010 August 12 (MJD 
55397 — MJD 55420). From these observations with a time span of ~ 23 days, we find 
that arrival times can be fit to a quadratic. From this fit, we calculate a spin frequency and 
spin frequency derivative of 1.8734(8) x 10~ 3 Hz and 6.42(6) x 10~ 12 Hz/s respectively. The 
residuals of the quadratic is further found to fit a circular orbital model with - sin % = 65 ± 3 
lt.s and an orbital period of 15.3 ± 0.2 days. 

Using this ^sini value, we find the mass function pf ) to be about 1.3M . An 

orbital 

orbital period of 15.3 days and the a spin period of 533.76s makes the source to be in line 
with the accretion powered pulsars with supergiant companions (Drave et al. 2012, Corbet 
1984). On the other hand, small mass function obtained from the circular orbital model 
should be an indication of a small orbital inclination angle. If the circular orbital model is 
preferred as the model fitting the residuals of the quadratic fit, this indicates that we observe 
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Figure 7. Pulse phase resolved variations of the spectral parameters with 5 phase bins for the first three RXTE— PCA 
observations. The pulse profile is plotted in 10 phase bins at the top panel. From top to bottom, we show the pulse profile, 
2—10 kcV unabsorbed flux, the photon index, the folding energy of exponential rolloff and reduced \ 2 of the fits. All uncertainties 
are calculated at the 90 per cent confidence level. For clarity, the data points are repeated for a cycle. 

the binary system nearly edge-on. Despite this configuration of the binary orbital plane, we 
do not observe a significant change in the Hydrogen column density which is not plausible. 

Alternatively the residuals of the quadratic fit can be explained by a torque noise strength 
of 6.8 x 10~ 18 Hz sec -2 . This value is quite consistent with other accreting X-ray binaries 
(Baykal & Ogelman 1993, Bildsten et al. 1997). Future observations is needed to understand 
the exact nature of the source. 

Initially, a double peaked pulse profile is observed in the light curves of the source. Then 
we find that one peak loses its intensity starting from the middle of the observation on 
MJD 55406.7 as the source flux continues to decrease after the burst. To study the energy 
dependence and temporal variability of pulse profiles, we construct pulse profiles with five 
different energy bands shown in Fig. |H We observe stronger pulses in the 8—13 keV and 
13—18 keV energy bands but generally the pulse shape shows no strong energy dependence. 
Double to single peak transition seen in this source might be due to a sharp decline in 
the intensity of the radiation coming from the fan beam, since the formation of fan beam 
strongly depends on the luminosity of the source as for EXO 2030+375 (Parmar et al. 1989) 
and GX 1+4 (Paul et al. 1997). However lack of significant energy dependence of the pulse 
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profiles makes the fan beam explanation implausible since fan beams are expected to be 
spectrally harder than pencil beams. 

In order to have a basic understanding of spectral variability of the whole observations, we 
study hardness ratios of the source. From the hardness ratio plots (see Fig. [5]), we suggest that 
the emitted radiation becomes softer as the source flux decreases during the ongoing decline 
of the outburst. The similar spectral softening with decreasing flux was reported before for 
1A 0535+262, A 1118-616, SWIFT J1626.6-5156, XTE J0658-073 and GRO J1008-57 using 
hardness-intensity relations (Reig & Nespoli, 2013). 

To extend the preliminary spectral analysis of Markwardt et al. (2010b), we also construct 
a single spectrum from all the available RXTE and Swift observations (see Fig. [6] and Table 
[2]). We find that adding an Fe line complex feature with a peak at 6.51 keV slightly improves 
the spectral fit. We discuss that this Fe line feature is more likely originated from the source 
as the galactic ridge emission is too weak to explain this emission alone. 

We also perform pulse phase spectral analysis of the source using the first three RXTE— PCA 
observations. From this analysis, we find that photon index and folding energy of the high 
energy cut-off vary with varying pulse phase (see Fig. |7]) . 
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